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Abstract 



First results of lattice QCD simulation on the nucleon tensor-charge 5q 
are presented. From the quenched QCD simulations with the Wilson quark 
action at (3 = 5.7 on a 16 3 x 20 lattice and on a 12 3 x 20 lattice, we find (i) 
the connected part 5q con , is determined with small statistical error, (ii) the 
disconnected part Sqdi S ., which violates the OZI rule, is consistent with zero, 
and (hi) the flavor-singlet tensor charge <5£(/x 2 = 2GeV 2 ) = 5u + 5d + 5s takes 
0.562(88) at f3 = 5.7, which is in contrast with the flavor-singlet axial charge 
AS = 0.1 - 0.3. 

PACS numbers: 13.88.+e, 12.38.Gc, 14.20.Dh 
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The parton structure of the nucleon at the twist two level is known to be characterized by 
three structure functions f±(x, /i), gi(x, /i) and h±(x, /i) with x being the Bjorken variable and 
\i being the renormalization scale (see, ]TJ). f± and g%, which represent the quark-momentum 
distribution and the quark-spin distribution in the nucleon respectively, can be measured 
by the deep inelastic lepton-hadron scattering. On the other hand, hi(x), which represents 
the quark-transversity distribution, could be measured in the polarized Drell-Yan processes, 
since it is related to the matrix element of the chiral-odd quark operator. Such experiment 
is planned using the relativistic heavy ion collider (RHIC) at BNL. Therefore, theoretical 
prediction of h\ is of great importance. Also, whether there is a large OZI violation in the 
first moment of h\[x) has particular interest since the first moment of g\[x) and that of 
the twist three structure function e(x) have been argued to have large OZI violation from 
experimental and/or theoretical point of view @^JQ]. In this paper, we will report first 
results on the first moment of h\(x) in lattice QCD simulations 0. 

The tensor charge of the nucleon is defined as 5q(fi) = / 1 [/i 1 (x, /i) — hi(x, fi)]dx, where 
hi(x) and hi(x) are the quark and anti-quark distributions respectively. One can also write 
matrix element of the tensor operator 

(ps | qia^q \ ps) = 2(s^p v - s„p^) 5q, (1) 

where p^ is the nucleon's four momentum, and s M is the nucleon's covariant spin- vector nor- 
malized as s 2 = — 1 p[. In the nucleon light-cone frame, 8q is interpreted as the "transver- 



sity" of quarks in a transversely polarized nucleon Q. On the other hand, in the nucleon 

q | ps). This implies that Sq 



rest frame with /x = and v = 3, 5q ~ (ps \ q 

\ a 3j 

gives an independent information on the spin structure of the nucleon from the axial charge 



Aq ~ (ps | g7 3 7 5 g | ps) ~ (ps | q 



'a, o ^ 



q | ps). 



y -o- 3 1 

As for Sq, approximate but model- independent inequalities such as | 5u \< 3/2 and 
I 5d \< 1/3 are known J7[. The non- relativistic quark model predicts 6u = Au = 4/3 and 
5d = Ad = —1/3, while the relativistic quark wave functions with non- vanishing lower 



components lead to 8u = 1.17 and Sd = —0.29 together with the inequality | 5q \>\ Aq 
There also exist estimates of 5q using QCD sum rules |||9j and a chiral quark model [[10 
but the uncertainties of the results are either large or uncontrollable. 

In lattice QCD simulations, one can treat both the connected (OZI preserviong) and 
disconnected (OZI violating) contributions even in the quenched approximation as has been 
demonstrated in [PLlT,pL2| - Also, one can estimate systematic errors due to the lattice approx- 
imation by combining simulations with different lattice spacing and/or lattice size. 



On the lattice, the matrix element of the tensor operator 0(t,x) = q 
extracted as 

E s= ±s(N s (t) £ t¥0)X O(^x)# s (0)> 



V a3 / 



q is 



R(t) 



Es=±(N s (t)N s (0)) 
const. + Z^Sq t (large t), (2) 



where Zt is the lattice renormalization factor for the tensor operator and N(t) is the nucleon 
operator projected to zero momentum. 

Following the works in ref. |]TT[ where similar simulations have been done, we have taken 
into account the points (i)~(iv) below. 

(i) We adopt the Wilson quark action and take N = {q t C~ 1 'y^q)q as a nucleon interpolating 
field. To enhance the nucleon signal, we use a wall source at initial time slice t = with the 
Coulomb gauge fixing at t = 0. Other time slices are not gauge- fixed. 

(ii) To avoid the contribution of the negative-parity nucleon propagating backward in time, 
the Dirichlet boundary condition for the quark propagator in the temporal direction is 
adopted. 

(iii) The connected amplitude is calculated by the conventional source method |13|]. To ob- 
tain the disconnected amplitude, the quark propagator with space-time unit source without 
gauge fixing is adopted. By this procedure, only the gauge invariant closed-loop is automat- 



ically selected after the average over the gauge configurations Note that t' = must 



be excluded in the sum in eq.(0) to make the procedure consistent. 



(iv) To compare the matrix element on the lattice with that in the MS scheme, we employ 
the tadpole- improved renormalization factor Zt(jicl). At fi = l/a with a being the lattice 
spacing, Z T in the Wilson quark action reads [|i~4f 



where K is the hopping parameter. The strong coupling constant ajjg(l/a) is evolved down 
by the 2- loop renormalization group from aj^g(ir/a) given 

by «is( ff / a ) _1 = Pa o l +0.30928. 



We carry out the quenched QCD simulation with the Wilson quark action at (3 = 5.7 
on a 16 3 x 20 lattice. In order to estimate the finite size effect, simulations at (3 = 5.7 on 
a 12 3 x 20 lattice are also performed. Simulation parameters and some of the measured 
quantities are summarized in Table |. Gauge configurations, which are generated with a 
single plaquette action separated by 1000 pseudo heat bath sweeps, are analyzed for three 
hopping parameters. The u and d quarks are assumed to be degenerate with K u = K d = K q , 
while the strange quark is assigned a different hopping parameter K s . The statistical errors 
of the data points are estimated by the single elimination jackknife procedure, and final fit 
of the hadron masses and R{t) are done by the x 2 -fitting- The lattice spacing a is extracted 
from m p a in the chiral limit. 

The connected and disconnected parts of R(t) at (3 = 5.7 on a 16 3 x20 lattice are shown in 
Fig. [I] at K = 0.164. Clear linear behaviors are seen for the u and d connected contributions 
5q CO n. in the interval up to t ~ 11 even for the smallest quark mass K = 0.1665. For t > 11, 
errors grow rapidly. The disconnected contribution 5qdi S . are concentrating around zero with 
errors comparable to the signal below t = 11. Since the disconnected part of R(t) does not 
show clear signal of a linear slope and we find no other sensible criterion, we take the same 
interval to extract Sqdis. in our linear fit. We adopt the same interval also for the fit of 
hadron masses. The linear fit of the connected part, the disconnected part and their sum 
are done separately. Our data on Sq is given in Table [TI]. 

For u and d quarks the fitted values are linearly extrapolated to the chiral limit m q a = 




(3) 



Here P is the plaquette expectation value and ckq is the bare lattice coupling constant |TJJ . 
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[\jK q — l/K c )/2 = 0. For the strange disconnected contribution, we first make linear 
interpolation to the physical strange quark mass m s a = (1/K S — l/K c )/2 at each fixed 
value of K q for nucleon, and then the results are again linearly extrapolated to the chiral 
limit m q a = {1/K q — l/K c )/2 = 0. The physical strange quark mass m s a is estimated 
by generalizing the relation m 2 a 2 = Am q a, obtained from the hadron mass results, to 
m 2 K a 2 = A(m q a + m s a) /2 and using the experimental ratio mx/m p =0.64. The errors of the 
sum of u, d and s contributions are estimated by quadrature. 

In Fig. H Su con ,, 5d con , and 5udi S . = 5ddi S . at (3 = 5.7 on a 16 3 x 20 lattice are shown as 
a function of m q a = (1/K g — l/K c )/2. The values linearly extrapolated down to the chiral 
limit are also shown with error bars. 

Our final results extrapolated to the chiral limit at (3 = 5.7 on a 16 3 x 20 lattice are 
summarized in Table The results for L = 12 (/3 = 5.7), which are done to estimate 
the finite size effect, are also shown in the table. To compare our data on 5q with that of 
Aq, results of a previous simulation for Aq in ref. [|TTJ at (3 = 5.7 on a 16 3 x 20 lattice are 
presented. 

From Table || and [ITJ, one finds the following features. 

1. Both the connected and total tensor charges at (3 = 5.7 on a 16 3 x 20 lattice satisfy the 
inequalities | 5u \>\ Au \ and | 5d \<\ Ad \ . This is different from the prediction of naive 
quark models | 5q \>\ Aq | mentioned before. 

2. We did not see clear linear slope for the disconnected parts and we could not make 
definite conclusion of its precise value from our simulation. Nevertheless, rather small error 
bars (~ 0.05) for 5qdi S . in our main simulation (L = 16 at (3 = 5.7) indicate that the OZI 
violation in 5q is small with a conservative upper bound | 5qd is . \< 0.1. Also, one notices 
that the disconnected part suffers a large finite volume effect. This is seen by comparing the 
data with smaller lattice sizes (L=12 at (3 = 5.7). The smallness of Sqdis. could be related to 



the C (charge conjugation)-odd and chiral-odd nature of the tensor operator qa liu 'y 5 q |L6 



3. Due to the small disconnected contributions, the flavor singlet tensor charge 
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5u + Sd + 5s is not much suppressed from its quark model value. This is in contrast to the 
flavor singlet axial charge AS which is known to be largely suppressed |§^|]. 

We have also done the simulation at (3 = 6.0 on a 16 3 x 20 lattice (a -1 = 2.271(9)GeV, 
La = 1.39 fm) and obtained qualitatively the same results as above. In the future, sim- 
ulations on a larger volume at (3 — 6.0 should be done to extract a definite value of the 
connected and disconnected tensor charges in the continuum limit as well as to reduce the 
finite size errors. 
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FIGURES 

FIG. 1. R(t) as a function of t at K = 0.164 and (3 = 5.7 on a 16 3 x 20 lattice. 

FIG. 2. Quark mass dependence of the tensor charges at (3 = 5.7 on a 16 3 x 20 lattice. Crosses 
with error bars denote the values extrapolated to the chiral limit. 
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TABLES 






5.7 




L 3 x T 


16 3 x 20 


12 3 x 20 


# of conf. 


1053 


306 


K 


0.16, 0.164, 0.1665 


0.16, 0.164, 0.1665 


fit range 


6 < t < 11 


5 < t < 10 


K c 


0.16925(3) 


0.16928(7) 


a' 1 (GeV) 


1.418(9) 


1.452(25) 


La(fm) 


2.23 


1.63 




0.2158 


0.2158 


m p a 


0.5429(36) 


0.5303(92) 


m N a 


0.8024(65) 


0.7985(144) 


m s a 


0.0896(12) 


0.0857(30) 


A 


2.752(13) 


2.746(33) 



TABLE I. Summary of simulation parameters. 
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If 
K q 


fill-con. 


A/7 

""■con. 


fiUdis. = Mdis. 




OSdis. 




= 5.7, 


16 3 x 20 






if s =0.1665 


0.164 


0.160 


0.1665 


0.948(31] 


-0.196(15) 


-0.044(69) 


-0.044(69) 


-0.035(58) 


-0.009(51) 


U.104 


n on a ( 1 1 ^ 
U.9y4(il J 


-U.ZZU(DZ J 


-U.Uz l (Zo J 


fl (AAA f QQ\ 

-U.U4i(oo J 


-U.Uz / (zo J 


-U.Ulo(zO J 


n i fin 
U.loU 


l.U I Z(OU J 


-U.Zol(Zl J 


-U.UUo(lZ J 


n noo/ 1 fi"\ 
-U.UZZ(lt) J 


n m Q/i a\ 
-U.Ulo(14 J 


-U.UUo(lZ J 


p — 0.1, 


1 93 v on 






1^ — n 1 fifi^ 


n 1 fizL 


n 1 fin 

U. 1DU 


0.1665 


0.825(49) 


-0.227(23) 


0.081(84) 


0.081(85) 


0.067(74) 


0.053(66) 


0.164 


0.919(20) 


-0.242(9) 


0.034(38) 


0.035(44) 


0.034(38) 


0.026(34) 


0.160 


1 034(10) 


-0 257(4) 


0.014(17) 


0.016(22) 


018(19) 


014(17) 






TABLE II. 


Measured tensor charges at /x = 1/a 






tensor charge 


axial charge (ref. []TTJ ) 


P 




5.7 




P 




5.7 


size 




16 3 x 20 


12 3 x 20 


size 




16 3 x 20 


Su 




0.839(60) 


0.822(83) 


Au 




0.638(54) 






0.893(22) 


0.760(39) 






0.763(35) 


8d 




-0.231(55) 


-0.159(75) 


Ad 




-0.347(46) 


5d con 




-0.180(10) 


-0.220(17) 


Adcon. 




-0.226(17) 


5(u,d)dis. 


-0.054(54) 


0.076(71) 


A(u,d)dis. 


-0.119(44) 


SSdis. 




-0.046(34) 


0.071(46) 






-0.109(30) 






0.562(88) 


0.733(121) 


AS 




0.18(10) 



TABLE III. Tensor charges extrapolated to the chiral limit together with the axial charge 
obtained in ref. |jTl]| . The flavor singlet tensor charge is defined as <5£ = 5u + 5d + 5s. The matrix 
elements are evaluated at fi = 1/a. 
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